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Abstract: We investigate the stochastic behaviour of long wavelength modes of light 
spectator scalar fields during inflation. When starting from a classical field value, the 
probability distribution for the spectator both spreads out and moves towards an equilib- 
rium distribution. We study the timescales for a mixed quadratic and quartic potential. 
The timescale of equilibration depends on the parameters of the model, and can be surpris- 
ingly large, even much more than thousands of e-folds. These results imply that the initial 
conditions for spectator fields are not automatically erased during inflation. Applying the 
results to the curvaton model, we calculate the probability distribution of the curvature 
perturbation and discuss 'typical' Universes. 
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1 Introduction 

Scalar fields can play a significant role in the evolution of the very early universe. The 
inflaton itself is an obvious example, since it drives the large-scale dynamics of the entire 
Universe. Scalar fields feature in many models motivated by particle physics, and are often 
called moduli fields. Even if the scalar fields do not drive inflation, they can still play 
a significant role in various physical processes. We call this type of field spectators, and 
denote them as a. The curvaton [1] is a particular example of a spectator field, which 
after inflation gives rise to the observed curvature perturbation £ and has been discussed 
extensively in the literature [2]. Other examples are provided by the quintessence field [3], 
and by the MSSM flat direction fields, which can be considered as light spectators during 
inflation with interesting fluctuation properties [4]. There are also models of inflation, 
such as assisted inflation [5] and N-flation [6] in which multiple scalar fields collectively 
contribute to the inflationary expansion. 

The interesting question is: what happens to the spectators during inflation, and what 
are their possible values at the end of inflation? We focus on light spectators, which are 
subject to inflationary fluctuations. The initial field values of the spectators may either 
have been determined by processes preceding inflation or by a phase transition during 
inflation (see for example [7]). In all cases, during inflation the long- wavelength modes of 
the light spectators will be subject to a stochastic evolution that can be described by a 
Langevin equation and the ensuing Fokker-Planck equation, which yields the time evolution 
of the spectator probability distribution [8] (see also e.g. [9, 10] for some early references 
on the subject). This is the approach we adopt in order to study the dynamics of spectator 
fields during inflation. 
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Formal solutions for the behaviour of a light, self-interacting field in a de Sitter back- 
ground have been given by Starobinsky and Yokoyama [11]. There are also analytical ex- 
pression for the stationary equilibrium solutions, which describe the equilibrium achieved 
after a sufficiently long period of inflation [11]. In the equilibrium limit, all information 
about the initial field configuration has been wiped out. To describe the capability of in- 
flation to wipe out information from initial conditions, it is necessary to investigate what 
period is long enough to attain a distribution close to the equilibrium one. Thus, we as- 
sume that the spectator has a classical initial value <7o at the beginning of the final period 
of inflation; whether this value is attained before inflation or during inflation is irrelevant 
for our purposes. We then solve the time evolution of the distribution, and ask how many 
e-folds of inflation it takes to get close to the equilibrium distribution. 

We concentrate on the spectator field potential V = m?(j) 2 /2 + A(/> 4 /4. In the limit 
A — > 0, the distribution is Gaussian and analytically solvable. However, the general case 
requires a numerical solution of the Fokker-Planck equation, which we show in section 2. 
We discuss the evolution of the probability distributions as a function of the parameters 
of the potential and the number of e-folds N. We determine both the relaxation time, 
which measures the rate by which equilibrium is approached, and the decoherence time, 
which measures the spreading of the initial delta-peak like distribution 1 . We also discuss 
the evolution of the lowest re-point correlators and pay particular attention to the transient 
evolution that takes place before the system has time to decohere. An important conclusion 
is that the equilibrium behaviour may only be approached after a large number of e-folds, 
because the transient evolution can easily take hundreds or hundreds of thousands of e- 
folds, depending on the parameter values of the spectator potential. 

After discussing the general evolution of spectator fields during inflation, we then focus 
on how this translates into a probability distribution of the curvature perturbation £ in 
the cases where the spectator field is a curvaton (section 3). We conclude in section 4. 

2 Behaviour of spectator fields 
2.1 The Fokker-Planck equation 

Let us consider a light (m <C H*) real scalar spectator field o evolving in a de Sitter 
Universe. Inflation is sustained by some other physics, not necessarily a scalar inflaton. 
For simplicity, we assume that the Hubble rate remains constant throughout inflation. 
After integrating out the short-wavelength modes k <C H* of the spectator field (using 
an appropriate window function), the evolution of the long- wavelength modes a can be 
approximated by a Langevin equation of the form 

where the random Gaussian noise £(t) has the correlator 

(Z(t)t(t>))=S(t (2-2) 
1 This is not decoherence in the quantum sense. 
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Figure 1. Quadratic potential: evolution of the probability distribution starting from a Gaussian 
(green) converging to the equilibrium distribution (red). The blue dashed line is the potential. The 
parameters used are A = 0, m/H* = 0.2, do/ = 1, w(0)/H* = 0.1 which give N re i = 75 and 
N dec = 37.5. 



As a consequence, the evolution of the probability distribution of the spectator a can be 
shown to obey a Fokker-Planck equation, which reads [11] 

H*, N) = ^jf"{o)P{o, N) + ^L V '(a)P'(a, N) + ||^>, N), (2.3) 

where the dot is a derivative w.r.t the e-folds N, and the prime is a derivative w.r.t a. Given 
the initial probability distribution, solving (2.3) yields the distribution for all N. In what 
follows we will assume a fixed classical initial value cjq, essentially a delta-peak distribution 
~ 5(a — o"o), which for numerical purposes will be taken to be a narrow Gaussian. 
Let us further assume a potential given by 

V{a) = \m 2 a 2 + \\o\ (2.4) 

where A < 1. In the N = oo limit, the probability distribution will reach a stationary form, 
given by 

V 2A eXP {~^J Kl / 4 



m 4 ir 2 
3A 



where K n (x) is a Bessel function and TV is a normalization constant. The equilibrium 
distribution (2.5) is independent of the initial condition. The question then is: how long 
does it take for the distribution to reach near equilibrium, and what is the rate by which 
the initial delta-function spreads out? 
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2.2 Quadratic potential 

Let us first consider the quadratic potential, i.e. the limit A = 0. This simplest case is 
solved analytically using the Gaussian ansatz 

Using (2.3) with arbitrarily narrow initial width, we determine the average cr c (N) and the 
width w 2 (N): 

(a)(N) = a c (N) = a c (0) exp ("|^) , (2-7) 
(a 2 )(N) - al(N) = W \N) = - ( - W 2 (0)) exp (-%£n) , (2.8) 



87r 2 m 2 \87r 2 m 2 / \ 3i/ 2 

where <r c (0) = o"o is the initial central value of the distribution. 
From the exponents we read off two time-scales (or iV-scales) 

= 3 S = M?- (2 - 9 » 

which we name the relaxation and decoherence time-scales, respectively. The relaxation 
time measures the rate of approach to the average field value at equilibrium, whereas the 
decoherence time measures the rate by which the initial narrow distribution broadens and 
spreads out towards the equilibrium width. 

In Fig. 1 we show the evolution of the distribution starting from gq/H* = 1 and 
w(0)/H* = 0.1, with m/iJ* = 0.2, so that iV re i = 2iVd ec = 75. We see how the initial 
peaked distribution (green) remains a Gaussian throughout and approaches the asymptotic 
form (red). The blue dashed line is the potential in some arbitrary normalization. 

2.3 Quartic potential 

The solution for the pure quartic potential (m = 0) is more complicated compared to 
the solution for the pure quadratic potential. For the quadratic potential, the probability 
distribution is Gaussian at all times, with a mean and width that evolve in time. For the 
quartic case, this is no longer true. The initially narrow Gaussian distribution becomes 
distorted before settling down to its asymptotic form. 

For the quartic potential we need to solve the Fokker-Planck equation numerically, and 
then extract the time-scales analogous to iV ro i and iVd ec m the quadratic case. For numerical 
stability reasons, the equation has to be solved using an implicit discretization. Also, a 
numerical solution limits the range of A and m 2 directly available to us; in particular, 
the small A region remains inaccessible numerically. Nevertheless, the results allow us to 
confidently extrapolate to small couplings. 

We plot the evolution of the probability distribution in figure 2 at different values of 
N. The initial distribution is given by a narrow Gaussian with u;(0)/ii/* = 0.1, a^/H* = 
1 (green curve) where m = and A = 0.003125. The grey lines are spaced by 20 e- 
folds. The initial Gaussian is transforming into the equilibrium distribution, which is not 



- 4 - 



0.2 







1 1 1 ' 

1 N = 




i 


1 N = 20 
i / 




/ 
/ 

\ i 


\ N = 40 /, 




\ i 
\ i 
\\ ' 


\ N = 6() k 

\ Im 
\ ^jtm 




\\ / 
si ' 
^ ' 

{§» ; 


mmmW, / 1 




/ 

^» / 

1 

\ i i 



o/H !t 



Figure 2. Quartic potential: evolution of the probability distribution starting from a Gaussian 
(green) converging to the stationary distribution (red). The blue dashed line is the potential. The 
parameters used are m = 0, A = 0.003125, ao/H* = 1 and w(0)/H* = 0.1 and the black lines are 
spaced by 20 e-folds. 
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Figure 3. Left: The connected 1, 2, 3 and 4-point correlators for A = 0.03125, m = 0. Also 
shown is the noise-less, pure slow-roll evolution (magenta). Right: The deviation of the correlators 
from their equilibrium values. When the lines becomes straight, the correlators are exponentially 
approaching their equilibrium values. The transient behaviour at small N depends on the initial 
conditions. 



Gaussian in this case. Superficially, the evolution of the probability distribution function 
appears similar compared to the quadratic case of figure 1. However, there are interesting 
differences, as we now discuss. 

To quantify the time dependence of this solution, we plot the connected 1, 2, 3 and 4- 
point correlators in figure 3. The plot on the right hand side of figure 3 shows the difference 
between the correlators and their equilibrium values. The vertical axis is logarithmic and 
thus we can see that the n-point correlators are well approximated by exponential laws, 
since they appear linear in this logarithmic plot. To find the relaxation and decoherence 
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Figure 4. Quartic potential: The dependence of the relaxation, decoherence and other timescales 
on A -1 / 2 . Parameters used are ctq/H* = 1, m = and w(0)/-ff* = 0.1. Note that the lines Ndec 
and N4 lie exactly on top of each other, as do the N3 and N re i lines (within numerical errors). 



times related to this exponential behaviour, we make a fit using the ansatz 

N 

C n (N) = C n (oo) - (C^oo) - C n (0))e-^ . 



(2.10) 



Here C n is the connected n-point correlator. Analogously to the quadratic case, we denote 
N x = N Tel and N 2 = N dcc . 

For different values of the coupling A we perform such exponential fits, also establishing 
the range of N for which the exponential behaviour is realised. This gives us figure 4. We 
have chosen to rescale the abscissa to A" 1 / 2 , to illustrate the dependence of these timescales. 
Not only are they straight lines, but the odd order correlators have the same scaling. The 
even order correlators have another scaling, but one is twice the other, in complete analogy 
with the Gaussian case (equations 2.7, 2.8). The timescales are given by 



N, 



rel 



N, 



dec 



Na 



11.3 
5.65 



(2.11) 
(2.12) 



Although the evolution of the n-point correlators is exponential for much of the evo- 
lution (i.e. looks linear in figure 3), during the first few hundred e-folds the evolution is 
non-exponential. This initial non-exponential behaviour we called transient behaviour. 
Here the distribution develops a 3-point correlator as it rolls down the potential. To a rea- 
sonable accuracy, we have found that the exponential behaviour begins around N = N dec 
in the evolution. We have also verified that by the time the exponential regime has been 
reached, the time-scales iV re i and N dec are independent of the starting distribution. 

Although both the initial and the equilibrium distribution are symmetrical about the 
mean, the distribution can be quite skewed for the transition period, which can last hun- 
dreds, even tens of thousands of e-folds, depending on A. The smaller A, the longer the 
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Figure 5. Mixed quadratic and quartic potential: The dependence of the N re i on A -1 / 2 and m 2 . 
Parameters used are ao/H* = 1 and w(Q)/H* = 0.1. The upper solid line is m 2 = 0; m 2 increases 
in the lower lines. Note that these lines would asymptotically meet the points on the A = line, 
which should be placed at oo. 



transition period. Therefore we conclude that even in the relatively simple quartic case, the 
light spectators are not generally well described by their equilibrium distribution. Rather, 
it would appear that for inflationary scenarios where inflation lasts only for 0(100) e- 
folds, the distributions of light spectators are dominated by their value before inflation and 
subsequent transient behaviour. 

2.4 Mixed quartic and quadratic potential 

For the full potential (m 2 > 0, A > 0) we again solve the Fokker-Planck equation and find 
an exponential regime for the behaviour of (a) and (o~ 2 ) c . Again we find that to a good 
approximation (i.e. to the numerical precision and the degree to which the evolution is 
strictly exponential), iV re i = 2N^ ec . In Fig. 5, we again show iV re i with the abscissa A" 1 / 2 . 
We note that the dependence is no longer linear in A" 1 / 2 . The solid black line is the pure 
quartic m? = result, which is a straight line. Increasing the mass to m? = 0.05/16, 0.05/4 
and 0.05 gives the dashed lines. The solid (magenta) vertical line gives the pure quadratic 
result, and should be placed at 1/A 1 / 2 = oo. The dashed lines will asymptotically approach 
the points on this A = line, for small enough A. When A decreases, the relaxation time 
becomes more and more dominated by the mass term. For smaller m 2 a similar behaviour 
can be found, where the relaxation time depends on A down to a certain value of the 
coupling, below which it approaches the A = limit at that particular m 2 . Thus, for either 
large m? or large A, the relaxation time is short. 

We emphasise that because of the noise term in the Fokker-Planck equation (the last 
term in 2.3), the relaxation time is not simply a function of the ratio m 2 /A as for the 
classical dynamics. 
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3 Initial conditions for the curvaton 



3.1 Curvaton probability distribution 

During inflation all light spectators obtain an isocurvature perturbation; the curvaton is 
a specific type of spectator which survives inflation and whose isocurvature perturbation 
is eventually converted to the observed adiabatic curvature perturbation £. Once inflation 
has ended and the inflaton has decayed, the curvaton starts to oscillate in its potential. We 
assume that the Universe is radiation dominated at that point. The relative amplitude of 
the curvaton field perturbation is initially negligible, but is then enhanced so that at the 
time of curvaton decay, the perturbation can be imprinted on the decay products as the 
dominant adiabatic perturbation. At the time of the decay, the curvaton energy fraction 
may either dominate over the inflaton decay products or be subdominant; this has an 
effect on the magnitude of the generated curvature perturbation. Here we will denote the 
effective decay rate of the curvaton by T. 

The simplest curvaton potential is a quadratic one. The perturbation amplitude and 
especially the non-Gaussianities of the perturbation may depend on the form of the poten- 
tial, as discussed in [12] . Here we do not aim at a detailed survey of the model parameters 
so that it suffices to consider the simplest quadratic case with A = 0. In that case, the 
curvature perturbation £ is given by 



if* 3p a 



' 37TCJ* 3p a + 4/9 



rest 



dec 3™ 



(3.1) 



where c* is the curvaton field value when the observable scales exit the horizon during 
inflation. The energy densities are evaluated at the time of curvaton decay and are given 

by 

* - 2W (f f ^ (5)) 3 ' ft - = M {wrf • (3 - 2) 

because they have the equation of state of matter and radiation, respectively. The curvaton 
slow-rolls until H(t) ~ m, which explains the factor oc (-^f) 3 ^ 2 ; Prcst is assumed to be 
responsible the Hubble expansion H* prior to <7* . We evaluate these quantities at curvaton 
decay, taken to be when 1/T ~ 1/H(t), where T is the effective decay width and radiation 
domination is assumed. This means that 

Pa\dec ~ m 1 ' 2 °l T 3/2 , PrestUc = 3M^T 2 (3.3) 

and thus we have 

C = — ^ (3 4) 

3 ^ + 4M| l( £) 1/2 

Inverting this, we find 
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Figure 6. The probability of ( < 10 -5 as a function of to, in the equilibrium distribution. If we 
require our Universe to have P(£ < 10~ 5 ) ~ 0.5, then the figure clearly shows a preferred range for 
m/H*. 



For fixed curvature perturbation f, there are two solutions for a*, denoted here as cr+ 
and cr~ , which correspond respectively to the dominant curvaton with r^ cc ~ 1 and the 
subdominant curvaton with r^cc <C 1. Note also that there are no solutions for a particular 
C unless 



1447r 2 Mg i r 1 / 2 C 2 

#2 m l/2 



< 1. 



(3.6) 



As an example, we now set T = 10" 15 GeV and H* = 10 10 GeV (M P i = 2.435 x 10 18 GeV 
is the reduced Planck mass). For these parameters, we need m > 70 TeV to get £ > 10 . 
The probability distribution of £ is 



P(C) = P[a-(C)] 




(3.7) 



where X 



1447r 2 M 2 / p 

— ff 2 \J mi anc ^ ^he ec i ua tions are valid only when 3.6 is satisfied. The 
distribution has non-trivial behaviour, given by an interplay between the condition 3.6, the 
width of P(cr) and the shape of C( c7 )- 



3.2 Typical Universes 

Observations indicate that our Universe has (£ 2 ) 



^10 5 ) 2 . Any curvaton model with 



£ > 10 is therefore ruled out. Models with £ <C 10 are not ruled out, but require 




1.5e+07 



Figure 7. Left: The evolution of probability distribution for £ for m/H* = 0.001 using iJ* = 10 10 
GeV, r = 10~ 15 GeV. Right: The probability of Q < 10~ 5 as a function of N for the same parameters 
as the left panel. 



an additional source of £. Figure 6 shows the probability of £ < 10 5 . If our Universe 
is 'typical', we would expect it to reside in the middle of the distribution, so that P(( < 
1(T 5 ) ~ P(C > 1(T 5 ) ~ 0.5. For small m, P(( < 10~ 5 ) = 1 because the condition (3.6) is 
not satisfied for £ > 10~ 5 . In this case, small values of m lead to r^ec <C 1 and small £. 
For large m, the distribution of P(cr) is narrow, and the a~ contribution dominates P (Q. 
This favours small C an d again gives P(£ < 10 -5 ) ~ 1. However, for intermediate masses, 
both the er" and the er+ contributions are important, and larger values of £ are probable. 
For these intermediate values of m, the width of P(cr) is such that the maximum of C( <J ) 
is favoured. Either larger values of H* or smaller values of V cause £ to increase, thus 
changing the range of masses that would give a 'typical' Universe. 

Let us now consider how the probability distribution of £ evolves as a function of the 
number of e- folds. We start with a very narrow Gaussian w 2 (0) — > around do = 0, 
and then let it evolve under the Fokker-Planck dynamics. The initially small width w 2 (N) 
increases towards the asymptotic value, which means that the a* and thus the C distribution 
broadens (figure 7 (left)). We used m/H* = 0.001, so iVdcc = 1.5 x 10 6 . We see the 
£ distribution widens in time, and moves to the asymptotic in a few million e-folds. In 
figure 7 (right) we show the integrated probability for ( < 10~ 5 , where the decoherence 
timescale is also obvious. 

4 Conclusions 

We have discussed the evolution of a spectator field in a de Sitter background, going beyond 
the formal solutions of [11]. We have given an explicit solution for the quadratic case, and 
solved the equation numerically for Gaussian initial conditions for a mixed quadratic and 
quartic potential. This is important because scalar fields can play a significant role in 
the evolution of the very early universe. Examples include inflaton fields, moduli fields, 
quintessence fields, MSSM flat directions and curvaton fields. 
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Evolution in the case of a quadratic potential is simple as one would expect: if the 
initial distribution is Gaussian, it will stay Gaussian for the remainder of its evolution. 
Furthermore, the mean and variance evolve according to exponential laws, and thus all 
information about the dynamics is encoded in two numbers, iV re i and A^ec- 

However, when one switches on interactions, the dynamics of the system become much 
more complicated. Even if the distribution starts out Gaussian, it will not stay that way. 
In fact, all higher order n-point functions are non-zero. On sufficiently long timescales the 
evolution of these correlators is exponential, and we have found the timescales that describe 
their evolution. Interestingly enough, reaching the exponential behaviour takes about iV"<j ec , 
before which the distribution evolves in somewhat more complicated ways. This transient 
behaviour can be surprisingly long, possibly many thousands of e-folds. Therefore, unless 
inflation lasts a very long time, this means that one really needs to investigate the full 
problem numerically if one wants to study spectator fields. 

Furthermore, we should point out that contrary to the general misconception that 
inflation erases information very quickly, the pre-inflationary conditions of spectator fields 
are erased very slowly, in many cases taking more than several thousands of e-folds. 

We have applied these results to the particular case where the spectator field is a 
curvaton and contributes to the curvature perturbation We calculated the probability 
distribution of £ for various masses (quadratic potential) in the equilibrium limit. We also 
showed how the probability distribution of £ evolves with time. These arguments could be 
used to quantify 'probable' Universes. It would be interesting to scan the whole parameter 
space of {m, T, H*, do, N} and discuss the likelihood of C = 10~ 5 . 
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